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SUMMARY 


The laminar compressible "boundary layer in the two-dimensional and 
axisymmetric stagnation regions has been analyzed to show the effects of 
the injection of a radiation absorbing foreign gas on an incident radia- 
tion field, and on the enthalpy profiles across the boundary layer. Total 
heat transfer to the stagnation region is evaluated for numerous cases 
and the results are compared with the no shielding case. Required absorp- 
tion properties of the foreign gas are determined and compared with 
properties of known gases . 


INTRODUCTION 


In some regimes of high-speed flighty the shock layer on a blunt 
body emits thermal radiation which is incident on the vehicle. This 
phenomenon becomes more prominent with increasing flight speed, diminish- 
ing altitude, and increasing nose radius. Kivel (ref. l) analyzes the 
radiation problem in the inviscid shock layer. His results show, for 
example, that the radiation heat transferred to the stagnation region 
exceeds the convective heat transferred if the nose radius of curvature 
is greater than about 1.7 feet and the vehicle is traveling at escape 
speed at 200,000 feet altitude. The surface is assumed to be black and 
to accept all of the incident radiation. 

The problem of shielding the vehicle surface from this thermal 
radiation becomes important . It is well known (ref. 2) that ordinary 
convective heating in the stagnation region can be greatly diminished 
by the injection of gas into the boundary layer. The possibility exists 
that the radiative heat transfer to the vehicle surface can also be 
if the gas injected into the boundary layer is opaque to 
radiation. Of course, the absorption of radiation by the opaque gas 
raises its temperature and thus the temperature gradient of the boundary 
layer at the wall, and therefore increases the convective heat transferred 
to the vehicle. Thus the question arises, is a net saving of heat trans- 
fer achieved by injection of the opaque gas, and if so, what absorption 
properties need this gas have in order to be effective? The purpose of 
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this paper is to study the effects of an opaque gas injected into the 
boundary layer on the combination of radiative: and convective heat trans- 
fer in the stagnation region of bluff bodies traveling at hypersonic 
speeds. However, an exact analysis of the interaction of a radiation 
field with a mixture of air and a foreign gas in the compressible laminar 
boundary layer is a very difficult problem. First, because it requires 
the solution of a set of nonlinear, coupled, partial differential integral 
equations. Secondly, in the boundary layer where large temperature gra- 


dients exist, it is doubtful that Kirchoff T s .'Law can be used to relate 

emission and absorption properties of the mixture locally. For these 

reasons, an exploratory analysis has been made in which simplifying 

assumptions make the problem more tractable. It is expected that the 

results of the approximate analysis will reta _n the significant qualita- ^ 

tive aspects of the actual physical behavior. If these results show that 4 

gains are achieved by use of an opaque gas, one would be encouraged to q 

investigate the problem further, experimentally or theoretically. q 


Programming the problem for solution on the electronic computer was 
done by Mrs. Yvonne Sheaf fer. 
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G P 

C 
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f, F 

— ► 

F 

C 

h 

I 

j 

k 

K 

l 


slope of lines in figure 4 and equation (35) 

specific heat at constant pressure 

Chapman-Rubesin function, ■ 

b e H-e 

coefficient of diffusion of foreign species into air 
dimensionless stream functions 
radiation flux 

ratio of total enthalpy to enthalpy exterior to the boundary 
layer 

static enthalpy of the mixture 

radiation intensity (integrated ever wave length) 
total enthalpy, h + — 
thermal conductivity 

absorption coefficient, defined by equation (7) 
dimensionless radiation intensity defined by equation (2l) 
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Pr 

Lewis number. — 

' Sc 

exponent in equation (l) , zero for two-dimensional case and 
unity for axi symmetric case 

pressure 

Prandtl number, 

k 

"total heat -transfer rate (convective and radiative) 

radius of cross section of body of revolution 

dimensionless transformed independent variable parallel to 
body surface (eq. (9)) 

Schmidt number, — 

’ PD 

velocity component parallel to surface 
velocity component normal to surface 
mass fraction of foreign species^ ^1 

distance along body surface measured from stagnation point 
distance normal to body surface 

dimensionless absorption coefficient defined by equation (26) 

constant in velocity relationship (eq. (17)) 

dimensionless transformed independent variable normal to body 
surface (eq. (lO)) 

coefficient of viscosity 

gas density 

stream function 

JL 

J2 

Superscripts 


r t r 


derivative with respect to the independent variable concerned 



Subscripts 


flow exterior to the boundary layer 
quantities at the wall 

foreign absorbing gas in the boundary layer 
air 


ANALYSIS 


The physical model chosen for analysis is represented by sketch (a) . 
The region between the shock wave and the body is divided into a shock 
layer of hot radiating air and a boundary layer consisting of a mixture 
of air and a foreign absorbing gas being Injected through and normal to 
the porous body surface. 



Sketch (a) 



H H 4=-> 


5 


The following assumptions are made regarding the properties of these 
regions. Additional assumptions will appear and he discussed where they 
are needed in the analysis. 

1. The radiant energy being emitted from the hot shock layer is 
incident on the outer edge of the boundary layer in beams of radiation 
of integrated (over wavelength) intensity I e (ref. 3 , eq. ( 5 ) and foot- 
note of ref. b, p. 53 ) , all assumed to be traveling normal to the wall. 

2. The air within the boundary layer is assumed to be transparent 
to the incident radiation and its emitted radiation is neglected in 
comparison with the incident radiation flux. 

3. The injected foreign gas, which will absorb a portion of the 
incident radiation before it reaches the body surface, has absorption 
properties which are assumed to be independent of the wavelength of the 
radiation (i.e., it is a gray gas). 

The mixture of air and the injected foreign gas in the boundary 
layer is assumed to be chemically inert. 

5. The surface of the vehicle is assumed to be cold and black; thus 
it absorbs all incident radiation which reaches the surface and it emits 
no radiation. 

Before proceeding to the mathematics of the problem, it should be 
noted that the assumption of a black vehicle surface is a limiting case. 
Ideally, if the surface of the vehicle were purely a reflector, it would 
not accept the incident radiation, and a nonabsorbing gas could be injec- 
ted for the usual convective heating protection. Indeed, in that case an 
absorbing gas would be undesirable in that it would trap radiation energy 
in the boundary layer and probably increase the convective heat transferred 
to the reflective surface by raising the temperature gradient at the wall. 
However, the characteristics of the actual vehicle surface will be between 
those of a perfect reflector and of a black surface and at some degree of 
surface absorptivity it will become desirable to inject an absorbing 
foreign gas. To simplify the mathematical problem, the analysis which 
follows considers the limiting case of the black surface. 

The partial differential equations describing the laminar compressible 
boundary layer of a binary mixture of gases in the presence of a radiation 
field are statements of continuity, the momentum theorem, conservation of 
energy, and diffusion of foreign species and are expressed, respectively, 
as 


^-(pur Q n ) + A( pYro n) = o 


( 1 ) 
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'duA /ouA <3 ( c)u\ dp 

pu K'dx-J pv \by) by \ by. ) dx 


( 2 ) 


pu -T— + pv 
ox 


( jl lb + A 

y ( x aV u2/2 1 

\Cp byj dy 

\ Pry by _ 


+ aj L eD I (hl ■ h2 ' 


div F 


( 3 ) 


dw dw d / _ dw\ 

pu S + pv ^-^V. pD V :=0 


( 4 ) 


The boundary conditions are 
at y = 0 


v = v w , u = 0, j = j w , w = W w (5) 

at y -► oo 

u -* Ue, j ■* j e , W - C (6) 

The exponent n in equation (l) is zero for two-dimensional flow 
and is unity for axisymmetric flow over a bod;y of revolution. The third 
term on the right-hand side of equation (3) can be derived from reference 5 j> 
page > 4 - 57 , and arises from molecular diffusion of the two gases. It will be 
neglected by saying C p2 ~ C p (and therefore hi » h 2 ) in order to 
simplify the treatment of the problem. 

The last term in the energy equation (3) is the rate of gain of 
energy per unit volume due to the radiation flux. For the case where the 
element of gas volume is exchanging radiant energy (by absorption and 
eramission) with all other element gas volumes inside and outside the 
boundary layer, the radiation flux is expressed by integration over all 
space. Then the energy equation (3) is a partial differential integral 
equation which is exceedingly difficult to solve. The problem will be 
simplified by neglecting the radiant energy emitted by the absorbing gas 
in comparison with the radiant energy abosrbed . An energy balance is 
still maintained of course; enthalpy and other forms of energy are simply 
not diminished by eramission of radiation. This simplification is justi- 
fiable if (l) either the absorbing gas is a good absorber and poor emitter, 
or (2) if the rate of energy eramission from the foreign gas is small com- 
pared with the rate of energy absorption. The second instance above can 
be visualized under some circumstances. It can be expected a priori that 
a large mass fraction of the dense boundary- layer region (near the wall) 


<! -4" r-i 
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will be absorbing gas, and a small mass fraction of the much less dense 
region (near the shock layer) will be absorbing gas. Therefore, the bulk 
of the absorbing gas is in the dense region which is cold compared with 
the shock layer. If the shock layer and boundary layer are behaving like 
black bodies and because they both emit energy at a rate proportional to 
the fourth power of their respective temperatures, then it is clear that 
the energy emitted by the shock layer (and thus absorbed by the foreign 
gas) is large compared with energy emitted by the boundary layer, and the 
latter can be neglected. For shock-layer emission deviating from black 
body conditions, the argument is weakened but not necessarily invalidated. 
For purposes of this exploratory analysis, the radiant energy emitted by 
the absorbing gas will be neglected in favor of the energy absorbed, and 
the results apply whenever this situation exists. 

In order to express the last term in equation ( 3 ) we consider that 
the radiation traveling in the negative y direction (toward the wall) 

Is absorbed by the foreign gas such that the fraction of the local inte- 
grated intensity absorbed in traveling the small distance -Ay is propor- 
tional to the local density of the absorbing gas and to the path length; 
or in the form of reference 6 (pp. 5 and 24) , for any fixed " x, 



(7) 


where K can be defined as an absorption coefficient. Because we have 
assumed that there is no scattering and have neglected emission in the 
boundary layer, there is only a y component of the radiation flux and 
thus the last term in the energy equation ( 3 ) is simply 


div F = — - ( 8 ) 

oy 

Equations (l) , (2), ( 3 ), (4), and ( 7 ) are to be transformed from x 
and y as independent variables to s and r\ making use of the Levy 
transformation (ref* 7 ) ? a stream function, several definitions and 
assumptions and some exterior flow relationships as follows. 

The Levy transformation is 


s = / Pe u e^e r o 2ndx 

Jo 


(9) 


n 


T) = 


Ue£o_ 

s/TsC 


y 


p ay 


O 


( 10 ) 
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A stream function is defined so that 


8\[r n 

= P^o > 
oy 


chjr _ 

8x 


= P vr o 


n 


(ID 


and the continuity equation (l) is satisfied. The following quantities 
are defined 


g(n) = j/Je 


ph 

PeM-e 


= C (ref. 8) 


( 12 ) 

( 13 ) 


where C is assumed constant and 


f (n) = W^e 


(lO 


from which 


f(tl) = i|r// 2ic 


( 15 ) 


In the a xi symmetric blunt -body stagnation region it is assumed that 


r G = x ( 16 ) 

At the outer edge of the boundary layer, :he external velocity is 
described by 


Ue = (Bx 


( 17 ) 


and it is assumed that 




e^e 




stagnation 


(18) 
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A 
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yields 


Use of equations (l6) , ( 17 ) , and (l8) in equations ( 9 ) and (10) 


s - ^Pe^e x ^( n + l) 

' 2(n + 1) 


(19) 


I n + 1 
PCpe^e , 


y 


P dy 


( 20 ) 


Thus for two-dimensional flow and axisymmetric flows, s is proportional 
to x 2 and x 4 , respectively, while q is proportional to y modified by 
compressibility. A dimensionless-radiation flux is defined by 


Un) 


1 

jePel i e u e r o 



( 21 ) 


Formally transforming equations (2) , (3) , (4) , and ( 7 ) to the new 
independent variables s and r\ by means of the newly defined quantities 
and the assumptions results in the following set of ordinary differential 
equations (with constant Pr and Sc) 


f m + ff ,T = 




(22) 


f o* T + 

F > ' 


_ 1 _ 

Pr 


- I 




( 23 ) 


f¥’ 



tt 


0 


(2k) 


l' - aMl = 0 


(25) 


where 


a = 


K / P . eP . eC 

7 ( n + !)P 


(26) 


The corresponding boundary conditions (5) and (6) transform to 
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at T] = 0 


f = f w , f' = 0, g = 0, W « ¥ w 


(27) 


at T] 00 


f f -> 1, g -* 1, W 0, 1 - l e ( 28 ) 

where f w is proportional to the injection rate (as will be shown 
subsequently) and the third of conditions ( 27 ) comes from the fact that 
the wall is much colder than the flow at the outer edge of the boundary 
layer. The last of conditions (28) is the boundary condition on equa- 
tion ( 25 ) , indicating that the local radiation intensity approaches the 
incident radiation intensity near the outer edge of the boundary layer* 


The set of differential equations (22) through ( 25 ) would be more 
tractable if the right-hand sides of equations (22) and ( 23 ) were zero. 
Then the equations would contain functions of only one independent vari- 
able q and therefore similarity solutions could be obtained. This 
approach will be used and the right-hand side of equation (22) will be 
neglected by virtue of the qualitative physical argument of reference 9 
(based on the fact that the surface temperature is much lower than T e ) . 
Now, of course, equation (22) with the right-hand side equal to zero is 
the familiar Blasius equation (ref. 10) where q is related to the 
Blasius i, and f(rj) and its derivatives are :*elated to the Blasius F(|) 
and its derivatives by 


T1 = V2 | 
f(n) = p(6)//2 
f’ (n) = p'U)/2 

= F "(t)/2s[2 


(29) 


Equations ( 29 ) cause the boundary conditions or 1 ( 11 ) to be compatible 
with the boundary conditions on F(|) (ref. ll) . Thus a solution of equa- 
tion (22) with the right-hand side equal to zero can be obtained at once 
from reference 11* 

Next, the right-hand side of equation ( 23 ) will be neglected because 
u e ^ « j e for the stagnation region in hypersonic flow. 


I_> |_i 4=-> 
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The boundary conditions ( 27 ) and ( 28 ) on the diffusion equation are 
mathematically sufficient, but are not very useful for initiating the 
numerical integration at rj = 0. For this reason, it is necessary to 
choose initial values of Vly and so that (l) a mass balance on the 

air at the surface is satisfied for a given injection rate and (2) 

W ( 00) -* 0. The assumption that the air does not penetrate the vail through 
vhich the absorbing gas is injected, leads to the following mass balance 
on the air at the wall 


A 

4 

1 

1 


Pw^W + Pw v w(l “ W V ) - 0 

which when transformed to the independent variable T) becomes 


(30) 


W f v = Sc f w (l - W w ) 


(31) 


Formally integrating equation (24) twice, making use of equation (22) , 
the boundary conditions ( 27 ) and ( 28 ) , and equation ( 31 ) > leads to 


" f Sc r 00 Sc 1 

1 - l(rfc) /fwSc J (f'O 

Equation ( 32 ) gives the correct value of so that W(co) 0. Equa- 

tion ( 31 ) gives the corresponding W' v used to begin the numerical 
integration of equation (24) . 

In the radiation absorption law (eq. ( 25 )), the absorption coefficient 
a (and therefore K) is assumed to be constant although in principle K 
could be any given function of temperature without destroying similarity. 


Method of Solution 


Equations (22) through ( 25 ) (with the right-hand sides equal to zero) 
subject to boundary conditions ( 27 ) and ( 28 ) were integrated numerically 
by the Adams -Moulton (ref. 12, p. 200, eqs. 6.6.2) predictor corrector 
method using the IBM 704 electronic data processing machine. Briefly, 
the sequence of machine computation was as follows: the quantities f 


V; 


f ! wu Swj ^w> P t > Sc, and a are specified [f v and f T V are obtained 


from reference 11 by use of equations ( 29 )) • Equation (22) (with the 
right-hand side equal to zero) was integrated numerically as a convenience 

r oo, v Sc 

for the machine computation and simultaneously the integral (f n ) dr) 
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vac evaluated for use in equation (32). Initial values W w and W' w were 
obtained from equations (32) and (31), respectively. Two values of g' 
were aruitrarily chosen. Each value of g ! w was used separately to 
interrate equations ( 24 ), (2p) , and (23) simultaneously (working in that 
order) . The values of g(uo) resulting from the two integrations were used 
to interpolate linearly to give an improved value on g'^ so that 
g(.::o) l. The integration was repeated until she boundary conditions were 
satisfied. 


In all of the numerical examples, Schmidt number and Prandtl number 
were taken to be 0.72 (and implicitly Lewis number is unity) . 


DISCUSSION OF RESULTS 


Let our first objective be to determine the radiation intensity at 
the wall in terms of that incident on the outer edge of the boundary layer. 
This Information is obtained from the solution of equations ( 22 ) , ( 24 ) , 
and (25) . The ratio of the intensity of radiation at the wall to that 
incident on the outer edge of the boundary layer is plotted as a function 
of a for two values of f w in figure 1 . The symbol f v is the value 
of the stream function at the wall and can be shown to be directly propor- 
tional to the mass injection rate by virtue of equations (ll) , (lp) , (9), 
and (lO) , which are combined to yield 


Pw v w ” 


L w 


Pe u e r o r Ve fc. 


(33) 


An injection rate corresponding to f W = -I.23I /J 2 leads at once to 
laminar separation (ref. 11 and eqs. (29)). This corresponds to an injec- 
tion rate of about 0.023 pound per square foot second for a vehicle with 
a nose radius of 1.7 feet, flying at escape speed at 200,000 feet altitude. 
Therefore the injection rate shorn in the figure corresponding to 
i'w = -l/s /2 is substantial. The equation for the two curves shown in 
figure 1 is obtained by writing the formal solution of equation (2p) 


is = e JET ” (3k) 

^e 

The Integral f ^° W dp has the values of 2.73 £nd O.968 for f v = -l/j 2 
and -l/2v/2, respectively. 

It is seen that at a given value of the at sorption coefficient a 
for the higher injection rate (f v = -1/^2 ) , much more radiation is 
absorbed in the boundary layer before reaching the wall. It can be noted 
that if a is about 2 or greater, almost all of the incident radiation 
is absorbed before reaching the wall. 


<£, -4" r— | H 
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Next , it is interesting to examine in detail how the radiation 
inuen^ity diminishes from the outer edge of the boundary layer to the 
wall. Pro! lies of the radiation intensity across the boundary layer are 
shown in figure 2 for various values of the absorption coefficient and 
- 0l ^ le stronger o lowing case = -l/\/ 2 . For increasing values of a, 
the region in which most of the absorption takes place tends to move away 
from the wall (tj = 0) . 


The absorption oi radiant energy within the boundary layer raises 
the local enthalpy level. The influence of absorption on enthalpy can be 
shown from solutions of the energy equation (23) which yields profiles of 
total enthalpy across the boundary layer* Approximately 40 examples were 
computed; some typical results are shown in figure 3. The lower curve 
corresponds 00 either no incident radiation or no absorption. For either 
increasing incident radiation intensity and a given absorption coefficient 
or increasing absorption coefficient and a given incident radiation 
intensity, the total enthalpy level at a given location (away from the 
wall) increases. 

The dashed curve corresponds to a very large specified l e . It 
"overshoots" (i.e., g > l) indicating that the local total enthalpy (and 
temperature if it is taken to be roughly proportional to total enthalpy) 
in the boundary layer exceeds that exterior to the boundary layer. This 
situation is not compatible with the physical model being analyzed. In 
the first place, energy cannot be transferred from a lower to higher 
temperature without an input of work. Secondly, for higher values of 
-le? the overshoot would be larger and would penetrate farther into the 
boundary layer, violating the condition that the bulk of the absorbing 
gas is at a temperature considerably less than that of the exterior flow. 
Finally, the neglected reradiation would tend to smooth out and diminish 
the overshoot. There is probably a region around the onset of overshoot 
where the solutions are of diminishing value. However, for lack of other 
ciiteria, solutions are presented up to, but not into, the overshoot region. 
A different analysis would be necessary to study the behavior of the high 
-le conditions. 

It can . be seen in figure 3 that the enthalpy gradient at the wall 
increases with increasing cl and -Ze. This gradient is of major signifi- 
cance because it is proportional to the convective rate of heat transfer 
to the wall. It is obvious that the reduction of radiant heat transfer 
le suiting from absorption tends to be offset by an increase of convective 
heat transfer. Before evaluating this situation, it is advantageous to 
look at the influence of absorption coefficient and incident radiation 
intensity on the enthalpy gradient at the wall in figure b, Two features 
of this figure are of particular interest. First, for a given a and f w , 
g'w appears to be linear with Z e ; that is. 


Z e -0 


a l e 


& w — (g'w') 


( 35 ) 



Ik 


where a = a(a) is the slope of the straight line corresponding to a 
given a. Secondly , the curves for a = 1.0 and 1.5 appear to lie out of 
order, between the a = 0.1 and 0.5 curves rather than above them. That 
is, for a given incident radiation flux (-Z e ) , increasing the absorption 
coefficient (a) first increases the enthalpy gradient at the wall (g ! w ) , 
but for higher values of a, the wall enthalpy gradient tends to diminish 
toward the no absorption value. Physically, this decrease in wall enthalpy 
gradient for higher values of the absorption coefficient indicates that the 
incident radiation energy is being absorbed farther out in the boundary 
layer and influences the wall enthalpy gradient less. 

The single point for a = 0.3 shown in the figure was computed to see 
if it lies in order, which it does (between the a = 0.1 and 0-5 curves) . 

So far, it has been shown that in shielding the wall from radiant 
energy the convective heat transfer to the wall is increased. It now must 
be shown whether or not a net saving in heat transfer results from injec- 
tion of an absorbing gas. Injecting a nonabsorbing gas into the boundary 
layer will diminish the total heat-transfer rate to the wall by reducing 
the convective part of the heat-transfer rate.?, as for example, is dia- 
grammed by the upper curve in sketch (b) . The injection of an absorbing 


Total 

heat- 

transfer 

rate 



Sketch (b) 

gas might be expected to diminish the total heat-transfer rate even further 
(such as shown by lower curve in sketch (b) ) . The total heat-transfer rate 
with injection of an absorbing gas will be compared with that with injec- 
tion of a nonabsorbing gas at the same rate 1 A will be compared with B) . 
The dimensionless total heat-transfer rate to a wall through which an 
absorbing gas is injected is obtained by sum ming the convective and ra dia- 
tive heat -transfer rates and dividing by j e ^(n + l)c3(p e n e ) stagnatin' 
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Thus 


q with absorption 
J~( n + l)Cp(p e P e ) stagnation 


‘'V 


-g'v 

Pr 


(36) 


V 

4 - 

L 

L 


Similarly, the dimensionless total heat-transfer rate to a vail through 
which a nonabsorbing gas is injected at the same rate is 


q no absorption 
^e \/( n + l)Cp"(p e M. e ) stagnation 


( g v)ot =0 

+ 

Pr 


(37) 


Comparing the heat-transfer rate with absorption to that without absorption 
at a given injection rate, by use of equations ( 36 ) and ( 37 ) , yields 


q with absorption _ ^w^ r _ §'w 


q no absorption iePr _ (g^)^ 


(38) 


The comparison is shown in figure 5( a ) tor an injection rate corresponding 
to f v = -l//2. The figure shovs that a definite advantage is obtained 
by injecting an absorbing gas. The effectiveness of this method of shield- 
ing against excessive heating increases vith increasing absorption coeffi- 
cient and increasing incident radiation intensity. For a = 1.5 and 
l e = - 0 . 5 ; the total heat -transfer rate at the vail is diminished by 
approximately 2 / 3 . The region to the right of these curves corresponds 
to the enthalpy overshoot condition and, because the analysis is not valid 
in that region, the curves are not extended beyond the overshoot boundary. 


It is worth noting that a combination of equations ( 38 ) , ( 35 ) , and 
(34) yields an expression for the curves of figure 5 (a) (although the 
curves vere not obtained by this expression) . It Is 


q vith absorption 
q no absorption 


a(a) + Pr e - (g T 


v) i e =0 


Z e Pr - (g’) 


v'lpt* 0 


(39) 


Using for example the value 2.73 for the integral (see below eq. (3*0), 

0.72 for Pr, 0.0714 for (g 'wh e =o l rom figure 4, the value a = 0.160 

corresponding to a = 1 yields the simple relationship 



(4o) 


q 'with absorption _ O.287 Z e - 0.0992 
q no absorption Z e - 0.0992 

which can be used to describe the appropriate curve in figure 5. 

It is also of interest to compare the total heat-transfer rate during 
injection of an absorbing gas to the total heat-transfer rate without any 
injection (A vs C in sketch (b) ) . The dimensionless total heat-transfer 
rate at the wall without injection is 


q no injection 

7 — — . . “ £ i — — 

(n + l)C 0 (p e u e ) stagnation 
The ratio of equation (36) to (4l) is 

q with absorptio n _ Z w Pr - g’-^ 

q no injection - TlPr- - f,~' 'i 

L e rr US 'w-'f^o 


(41) 


( 42) 1 


The comparison is shown in figure 5 (b) . The injection rate for the absorp- 
tion condition corresponds to f w = -l// 2 . Here, of course, the effect is 
very pronounced. . Examination of the figure si lows first of all that it is 
advantageous to inject a gas whether it absorbs radiation or not. As seen 
before, the heat-transfer rate will be diminished more by injection of an 
absorbing gas. However, for low levels of incident radiation intensity, a 
large absorption coefficient is not much more effective than a small 
absorption coefficient in reducing the total heat-transfer rate. Again, 
for laige values of incident radiation intensity, large absorption coeffi- 
cients are very advantageous. In particular, for a = I.5, the heat- 
transfer rate is only 1/5 that for no injection at all for all values of 
incident radiation intensity, and less than 1/2 that for injecting a 
nonabsorbing gas if Z e = -O.5. 

Finally, something should be said about the actual absorption 
coefficient K. The required K is dependent on the flight condition 
by virtue of its relationship with a (eq. ( 26 )). For a flight speed of 
34,000 feet per second at 175,000 feet altitude, and body-nose radii of 
1 to 10 feet, the required K is of the order of 10 5 to 10 7 ft 2 /slug for 
values of a from 0.1 to 1 . 5 * Little is known of the absorption proper- 
ties of even ordinary gases. Measurements by Eckert and others (ref. 5, 
pp. 304-386) can be interpreted to give an approximate "grey" K for 
purposes of this discussion. It appears that carbon dioxide has an 
absorption coefficient of the order of 10 sq ft per slug at ordinary 
pressures, for temperatures up to l 600 ° C. The coefficient for water 

^Equations (4-2) , ( 35 ) > a nd (34) could be combined to give the heat- 
transier comparison of injection with absorption to the no injection case 
corresponding to equation (39) . 
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vapor is roughly an order of magnitude higher than that of carbon dioxide, 
and so the gas we want to blow into the boundary layer should be 3 to 5 
orders of magnitude better than water vapor in this regard. 

It is encouraging to find that at least one gas has a high absorption 
coefficient. That gas is cesium vapor, which has an absorption coefficient 
of the order of 10 s square feet per slug at 0° C and 1-mm Hg in the wave- 
length region 2000 to 3000 A (which can be shown from ref. 13 ) . Demon- 
strably then, such opaque gases exist, and the problem is to find a 
suitable one. 


CONCLUDING REMARKS 


There are several notable results of this exploratory analysis of 
the effects of the injection of an opaque gas in the stagnation region of 
a blunt body traveling at hypersonic speed. At a black vehicle surface 
the reduction of the radiation heat-transfer rate by radiation absorption 
in an injected opaque gas is accompanied by an increase in the convective 
heat- transfer rate. However, the net effect is that a saving in total 
heat-transfer rate (radiative plus convective) of as much as 2/3 can be 
achieved by injecting the absorbing gas into the boundary layer. The 
opaque gas must have a high absorption coefficient (3 to 5 orders of mag- 
nitude higher than water vapor) in order to effectively reduce the total 
heat-transfer rate to the body. It is pointed out that under some condi - 
Hons, cesium vapor has an absorption coefficient 6 orders of magnitude 
higher than water vapor . It would be worthwhile to search for gases that 
have , . for actual flight conditions of interest, suitably high absorption 
coefficients as well as desirable injection properties. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., May 17 , i 960 
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Figure 3**“ Influence of* radiation absorption on enthalpy profiles for 

*V = - 1/jl 
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Figure p.- The effect of gas injection at a rate corresponding to f v = - l//2 on the total heat 

transfer rate. 










O 05 
XI CM 

H M 

i 

c Q 

Q> < 

85? 

X 2 


U P 



2 2 


w c £ <u 

w r- ca -*j -r-t „ 

5 - rt **-• 

0) g S % C 

^ 5 3 5 © 

H X3 ^ P j 


< 

z 


1 


•'if 

o 

■*H 

- £ 

» >> 


1 


g & 

iMo 


•§ 


® s § 
ss s, 
**3 | 
« * s 
JS * 

O rn 5 
33 « ^ < 

^ • * 
Xl -tH W 

t" -cl a B 

h a> -3 § 

« > S i 

>9 ^ i <U 

<SS2 

S£ 0 ! 


o o> 

JS CM 

H ^ 


■§£ 
oT < 

83 

K 2 
~ P 


< co 

ohst 
i8«°s 

o O p 


on 


2g§ 


Q 5 




■d ^ 2 o 

< So MO 
S£z|S 

S-fcRlS 


r* ^ " 

S?§SB 


0> O 

' ' w 


?OS 3? 
Sgfi^c 

o n ^ co 

-9z "5 

sl|gl £ 

CO o w H • 

Q &*!" 1 “ 

“<o«ao 

< g 3 2 P «r 

W JH Q O C O 

Figgis. 




O 05 
CM 

H « 


oT <J 

852 

P 2 
~ P 


< 

< 

z 


^co 
OHS 

c 5 W CO 
o O p i 

SPa do 
2<«w 

S88 S 
1251* 

^ CQ 05 

< ^o 

<U ® ^ 

O P P £ = 

&§sgg 

ii| 

DO o ^ 

U 4 CO 

•llil* 

M C ffl K 3J 

co o 53 H _ , • 

Q <3 2 ^ ^ 

”<QO 0q 

< §22^ 
S^qo c 

£Vc SSI 



■ § 

.2 P 

&g S 

c ® ts 5 01 

g DO « ^ O 

22^1“ 

0 oi3 rt 5 

5 rt o « 2 

_, «0 ** 0> QJ 

§3 t 32 

g IH <H O B 

°2*«o 

DO ft d 3 C 

“ " 5 

&3ii5 

<U P A S3 P 

Lj C _j Si 

.2 *3 St. 

q> o m-h m 

1 c Si 

■S ° ^ g s 

S T3 S ^ O 

■S s « S; « 

C T3 >> GO S 

s 2 I C 2 

^Cd^oS 

rt O o o; <u 

i* p .a u u 





